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ABSTRACT 

Proving the magnetic configuration of solar spicules has hitherto been difficult due to the lack of spatial 
resolution and image stability during off-limb ground-based observations. We report spectropolarimetric ob¬ 
servations of spicules taken in the He I 1083 nm spectral region with the Tenerife Infrared Polarimeter II at the 
German Vacuum Tower Telescope of the Observatorio del Teide (Tenerife; Canary Islands; Spain). The data 
provide the variation with geometrical height of the Stokes I, Q, U, and V profiles whose encoded information 
allows the determination of the magnetic field vector by means of the HAZEL inversion code. The inferred 
results show that the average magnetic field strength at the base of solar spicules is about 80 gauss and then it 
decreases rapidly with height to about 30 gauss at a height of 3000 km above the visible solar surface. More¬ 
over, the magnetic field vector is close to vertical at the base of the chromosphere and has mid inclinations 
(about 50 degree) above 2 Mm height. 

Subject headings: Sun: chromosphere — Sun: filaments, prominences 


1. INTRODUCTION 

Spicules are thin and very dynamic needle-shaped struc¬ 
tures, best seen at the solar limb, whose magnetic proper¬ 
ties are not well constrained to date. A detailed study of 
spicules require very high spatial and temporal resolution 
observations, being the main reason behind their poor char¬ 
acterization. Modern solar space-based observatories, such 
as the Japanese satellite Hinode ( [Kosugi et al. 20071 and, 
more recently, NA SA’s Interface Reg ion Imaging Spectro- 
graph (IRIS; |De Pontieu et al.||2014| l have improved some 
of the spicules properties by usin g fast cadence imaging 
and spectroscopy, respectively (e.g., De Pontieu et al.1|2007 


|Pereira et al. 2012; Tavabi et al.|201 1 and references therein). 

However, their magnetic configuration is not yet well con¬ 
strained. 

The magnetic field in spicules has been determined using 
direct analyses of spectral lines in polarized light, mostly us¬ 
ing the He I 1083.0 nm triplet and the He I D 5 587.6 nm 
line, or indirectly by using wave properties (Zaqarashvili et 
|al.|2007| |Kim et ak|2 0081 Ve rth et al.|2011) . From the inter- 
pretation of spectropolarimetric observations, we know that at 
chromospheric heights, spicules harbor field strengths of the 
order of 10 gauss and inclinations of about 35° with respect to 
the local vertical, although stronger magnetic fields were also 
found in spicules ([Trujillo Bueno et al.|2005| Lopez Ariste & 


|Casini|2005| |Ramelli et al. 2006). Later, |Centeno et al.| fl?010) 

carried out a 45 minute time series in sit-and-stare mode in 
the He I 1083.0 nm triplet, with a spatial resolution of about 
2", and at two different heights above the visible solar sur¬ 
face. They interpreted the data with the HAZEL (from HAnle 
and Zeeman Light) code (Asensio Ramos et al.|2008) and ob¬ 
tained field strengths up to 50 gauss and field inclinations of 
35°-40°. 

In this Letter, we present observations in the He I 1030 nm 
triplet obtained under excellent seeing conditions and report 
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Halpha channel 



FIG. 1.— Image taken from the H a slit-jaw camera at the VTT. In the 
image it can be distinguished a sunspot (right side), where the adaptive optics 
system of the German VTT was locked, and the target of TIP-II, the spicules. 
The vertical line represents the spectrograph’s slit. In the image it can be seen 
that the spicules have a dominant orientation. 


the variations of the magnetic field vector of solar spicules 
with geometric height. We describe the observations in Sec¬ 
tion 2. In Section 3 we briefly explain the diagnostic tech¬ 
nique, and in Section 4, we present the results and the discus¬ 
sion. 


2. ANALYSIS OF THE OBSERVATIONS 

The observations were taken on 2013 Apr 23 with the TIP- 
II instrument ((Collados et al.||2007|) installed at the German 
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VTT of the Observatorio del Teide (Tenerife, Spain). Dur¬ 
ing the observations, the Kiepenheuer Adaptive Optics Sys¬ 
tem (KAOS, von der Luehe et a l.|2003| l of the VTT was con¬ 
tinuously locked on a sunspot located at the east limb, dur¬ 
ing good seeing observing conditions. The slit of the spec¬ 
trograph was situated crossing the solar limb with an angle 
of ~ 45° between the slit and limb. Then, an area of about 
9" (25 slit steps with 0"36 step size) was scanned. For each 
slit position, the TIP-II camera recorded the four Stokes pa¬ 
rameters around the 1083.0 nm spectral region with a spectral 
sampling of 1.1 pm/px. The scanned spectral region contains 
the chromospheric He I 1083.0 nm triplet, the photospheric 
Si I 1082.70 nm line, and an atmospheric water vapor line 
at 1083.21 nm. The length of the spectrograph slit is 80" 
and the spatial sampling along the slit is 0"17. The slit posi¬ 
tion of the first scanning step can be seen in the H Q slit-jaw 
image (Fig. ITJ). The total integration time per slit position 
was 10 seconds (10 accumulations and 250 ms exposure time 
per modulation state). We estimate that the spatial resolution 
could be in the 0."7 - 1" range. These data were subjected 
to the typical data reduction process that includes dark cur¬ 
rent, flat-held, fringe correction, polarimetric demodulation, 
and polarimetric cross-talk correction. 

Figure [2jshows an example of TIP-II full Stokes vector ob¬ 
servations in the He I triplet at 1083.0 nm. In particular, we 
represent a single slit data taken with the TIP-II camera. The 
first panel shows the variation of the Stokes I profiles of the 
He I triplet, which is seen in emission in off-limb observa¬ 
tions, with height above the limb. The slit was crossing the 
spicular material in an inclined angle with respect to the so¬ 
lar limb. Thus, the variations of the intensity with height in 
Stokes I correspond to spicules, extending up to 8 Mm. Inter¬ 
estingly, spicules can be clearly distinguished in the spectra as 
a result of to the excellent seeing conditions and high spatial 
resolution of the data. In the wavelength dimension, Stokes 
I shows the fine structure of the He I 1083.0 nm triplet: a 
weak blue component at 1082.9 nm and the stronger, two red 
blended components located at about 1083.03 nm. 

In Stokes Q and U, the signals are clearly seen above the 
noise level. These signals are due to atomic level polariza¬ 
tion, that is, population imbalances and quantum interference 
between the magnetic sub-levels of degenerate atomic levels 
generated by anisotropic radiation pumping and the Hanle ef¬ 
fect. In our data, the positive reference direction for Stokes 
Q is parallel to the nearest limb. Thus, the presence of sig¬ 
nal above the noise in the Stokes U panel readily indicates 
the presence of a magnetic field inclined with respect to the 
local vertical direction (which represents the axis of symme¬ 
try of the incident radiation) since, according to the Hanle 
effect theory, the presence of a magnetic field modifies the 
atomic level polarization. Interestingly, in the data of Stokes 
U we can see that there is a change in the sign of the profile 
(from white, positive, to black, negative) indicating that there 
are changes in the magnetic field azimuth with height in so¬ 
lar spicules. The physics of the Hanle and Zeeman effects in 
this triplet is described in the papers by Trujillo Bueno et al. 
(2002), Trujillo Bueno & Asensio Ramos (2007), and in the 
book of Landi DeglTnnocenti & Landolfi 1 2004] ). 

Finally, Fig. [2] (right panel) shows that there is a consid¬ 
erable amount of Stokes V signals in the data well above 
the noise level, i.e., circular polarization profiles dominated 
by the Zeeman effect, with two lobes of opposite sign. The 
mere presence of Stokes V signals allow us to determine the 
magnetic flux (through the Zeeman effect). Given that Stokes 


Q and U are dominated by atomic level polarization and the 
Hanle effect and that these signals allow us to determine the 
field orientation, we can determine the strength of the mag¬ 
netic field vector rather accurately from the combined inter¬ 
pretation of the Zeeman and Hanle signals. 

3. INVERSION OF THE STOKES PROFILES 

In order to interpret the intensity and polarization profiles 
that we measure in solar spicules we need to assume a suitable 
model. In our case, we assume a constant-property slab model 
with a given optical depth r (Trujillo Bueno et aL|[2005j). 
We assume that the spicular material of the slab, located at 
a given height above the limb, is equally illuminated by the 
unpolarized, limb-darkened photospheric radiation field pro¬ 
ducing atomic level polarization. We also assume a constant 
thermal width of the line (modeled with a line damping and a 
Doppler broadening parameters) and a deterministic magnetic 
field vector. For computing the Stokes profiles, we take into 
account the joint action of atomic level polarization and the 
Hanle and Zeeman effects in the incomplete Paschen-Back 
effect regime. The slab model takes into account atomic level 
polarization for all the involved atomic levels and therefore 
it naturally explains the selective absorption/emission mecha¬ 
nisms we detect in the profiles. In the slab model, the emer¬ 
gent Stokes vector is: 


I - e ~ K * T T 

x — e: A s , 


, + K" 


_1 (l-^ T ) 
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where S = e/?)i and K* = K/r/i, with e, rji, and K the emis¬ 
sion and absorption profiles for Stokes I and the propagation 
matrix, respectively. In off-limb observations, I sun = 0. To in¬ 
vert the observed profiles, we use the HAZEL code. The code 
recovers seven parameters from the inversion of the observed 
Stokes profiles: the optical depth of the line at the central 
wavelength of the red blended component r re( j, the line damp¬ 
ing, the Doppler width, the bulk velocity of the plasma, the 
magnetic strength B, inclination 9, and azimuth x of the mag¬ 
netic field vector. For reference, the inclination of the mag¬ 
netic field vector is given with respect to the solar vertical (90 
degrees means parallel to the solar surface), while the azimuth 
of the magnetic field vector has its origin in the line-of-sight 
direction and is measured counterclockwise (zero degree az¬ 
imuth and 90 degree inclination corresponds to a magnetic 


field vector pointing to the observer; see also Fig. 3 in Orozco 


Suarez et al. 2014). The height above the limb, which fixes 


the degree of anisotropy of the incident radiation field, is cal¬ 
culated from the data. A detailed description of the HAZEL 
code and the underlying phy sics can be found in the pap ers by 
Trujillo Bueno et al. (2002), Trujillo Bueno & Asensio Ramos 
(2007), and Asensio Ramos et al.|(|2008|). Finally, to improve 
the signal-to-noise ratio of the observations, the data were 
down-sampled spectrally and spatially along the slit direction, 
yielding a final spectral and spatial sampling of 4.4 pm and 
0"51, respectively. 


4. RESULTS AND DISCUSSION: VARIATION OF THE MAGNETIC 
FIELD VECTOR WITH HEIGHT 

An individual fit (solid lines) to prototypical profiles (dots) 
of the He I 1083.0 nm triplet is shown in Fig. [3] The Stokes 
I profile is normalized to unity and Stokes Q, U, and V to 
the peak value of Stokes I (before normalization). The fit is, 
overall, very good. In this particular case, the field strength 
we obtain from the HAZEL inversion is B = [54,51] gauss, 
the field inclination 9 b = [42,79]°, which corresponds to a 
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Fig. 2.— Stokes I, Q, U, and V maps representing a single slit position of the full TIP-II map. X-axis represents the wavelength and the Y-axis the height over 
the solar limb. The Stokes I panel shows the fine structure of the triplet (in emission), with a weak blue component at 1082.90911 nm ( 3 Sj— 3 Po) separated about 
0.12 nm from the other two components at about 1083.02946 nm ( 3 Si— 3 Pi and 3 Si~ 3 P 2 ), which are indistinguishable. The two components can also be seen in 
circular polarization. In Stokes Q and U, only the red component can be seen above the noise level although just above the limb a slight blue component signal 
in Stokes Q stands above the noise. This tiny blue component is due to the selective absorption mechanism. The positive reference direction for Stokes Q is the 
parallel to the solar limb. 


Stokes I 



Stokes Q 


Stokes U 




Stokes V 



Fig. 3.— Observed (dots) and best-fit (solid line) He I 1083.0 nm triplet Stokes profiles corresponding to a prototypical spicular profile. Stokes I is normalized 
to unity, while Stokes Q, U, and V are normalized to their Stokes I maximum peak value. Different color codes represent different solutions. The positive 
reference direction for Stokes Q is the parallel to the solar limb. 
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Fig. 4.— Variation of the magnetic field vector (strength, inclination, and azimuth, bottom panels) with height above the solar limb. The results correspond to 
the solution that is compatible with a magnetic field vector align with the solar spicules, i.e., the vertical solution. The error in the determination of the height 
above the limb is given by the spatial resolution along the slit. In our case, it is twice the resolution element, i.e., ±740 km. 


field vector pointing outward from the solar surface, and the 
field azimuth \'h = [1,41]°. Note that the inferred vector mag¬ 
netic field from the analysis of the He I 1083.0 nm triplet is 
subjected to the so-called 90°ambiguity of the Hanle effect 
and to the well-known 180° azimuth ambiguity (see Sect. 4.2 

in |Orozco Suarez et al.|20f4] |Martinez Gonzalez et al.|2015| 

and more references therein). The first one is associated with 
the Van Vleck angles (54.74°and 125.26°) and tell us that the 
magnetic field inclination may have two ambiguous solutions 
lying at both sides of any of the Van Vleck angles. The two 
solutions we quoted above represent the 90° ambiguous solu¬ 
tions. One of them is compatible with a magnetic field vector 
rather inclined with respect to the solar vertical and another 
with more vertical fields. In both cases, the inferred mag¬ 
netic field strength and its azimuth change. Regarding the 
180° azimuth ambiguity, in the case of 90° scattering geome¬ 
try, it implies that the solutions = 180° - 9b with Xb = ~ Xb 


produce the same polarization signals, so they are ambiguous. 
In practice, these solutions correspond to a vector magnetic 
field pointing towards the solar surface. Finally, the inferred 
optical thickness at the central wavelength of the red compo¬ 
nent does not change in both solutions and is, in this particular 
case, T re d = 1.3. 

In Fig. [4] we show the variation of the magnetic field 
strength, its inclination, and azimuth with height above the so¬ 
lar limb. Between the two 90°ambiguous solutions, we have 
chosen the one that is compatible with more vertical magnetic 
fields. The reason is twofold. First, in the more vertical so¬ 
lution, the inclination of the projection of the magnetic field 
vector onto the plane of the sky (the plane perpendicular to 
the line-of-sight direction) is in line with the average inclina¬ 
tions that can be measured in monochromatic observations of 
spicules, such as that of Fig. IT] Second, for a long time now, it 
has been assumed that spicules are plasma structures standing 



























4 


D. Orozco Suarez, et al. 


nearly vertical on the solar surface. Such assumption has been 
adopted for their MHD modeling. Hence, we consider that the 
vertical solution to be the one that makes sense in the current 
understanding of spicules (e.g.,[Becker s] 19681 . Nevertheless, 
note that from the point of view of the He I 1083.0 nm triplet 
modeling, the two solutions are equally probable. 

The inversion of the data shows that in solar spicules, the 
magnetic field strength falls rapidly with height, going from 
80 gauss close to the visible limb to about 30 gauss at 3 Mm 
(on average). Above 3 Mm the dispersion increases with an 
rms of about 30 gauss and a mean of 40 gauss. The high dis¬ 
persion might be associated with the noise in Stokes V. These 
values agree with previous measurements that were in the 
range of the tens of gauss for solar spicules (|Truj illo Bueno 
|et al.|2005]|Centeno et al.|2010| ). 

In the case of the field inclination, at the bottom of the atmo¬ 
sphere, it is close to 0° (vertical fields) and increases sharply 
to ~ 40° above 1.5 Mm. The rms for the inclination also 
increases, being about ~ 5° on average, with the rms being 
less significant than for the field strength. The inclinations we 
infer agree well with what it can be tentatively derived from 
monochromatic images. To compare with previous observa¬ 
tions, we need to calculate the inclination of the magnetic field 
projected onto the plane of the sky. The results are shown in 
Fig. 0 The mean inclination is about 20° with an rms of 11°. 
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Fig. 5.— Variation of the inclination of the projection of the magnetic field 
onto the plane perpendicular to the LOS (plane of the sky). The dots repre¬ 
sents the inferred values and the asterisks and bars represent the mean and 
the rms, respectively, for a 500 km bin size in the horizontal axis. 


Regarding the field azimuth, it changes from -50° to 50°, 
and the dispersion also increases with height, with an rms of 
~ 55° above 3 Mm. This means that the field vector is point¬ 
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ing toward the observer at almost all heights (azimuths from 
-90° to 90°). Very high in the atmosphere there are positions 
where the magnetic field points away from the observer. This 
behavior of the azimuth can be well understood if we have a 
look at the H„ slit-jaw image of Fig. IT] where it seems that 
the spicules contained in our data set nave a preferred orien¬ 
tation instead of random orientations. This may be due to the 
proximity of the active region. 

5. CONCLUSIONS 

In this Letter, we have presented novel spectropolarimetric 
observations of solar spicules taken in the He I 1083.0 nm 
triplet. We succeeded in partially seeing some of the fine 
structure of spicules, which is visible up to 8 Mm height. We 
have detected clear linear polarization signals associated with 
atomic level polarization and the Hanle effect, while Stokes V 
is clearly due to the Zeeman effect. The data have allowed us 
to determine the magnetic field vector variation with height. 
We find that the field strength in solar spicules decreases from 
80 gauss at the visible limb to 30 gauss at about 3000 km, 
with a large scatter at greater heights, due to noise and/or av¬ 
eraging along the LOS. No strong variations of the field ori¬ 
entation are found, although the field is clearly more vertical 
at low heights and then goes close to 40 degree at 1500 km. 
These results represent a step forward in the determination 
of magnetic fields in chromospheric spicules and should be 
taken into account in their numerical or theoretical modeling. 
It would be interesting to record similar observations at higher 
spa tial resolution, e.g., with the GREGOR telescope ( |Schmidt| 
|et al.|201~2] l, and compare with low resolution results and see 
if there is any difference in the magnetic field configuration 
inferred for spicules recorded near active regions, in the quiet 
Sun and in coronal holes ([Beckers 1968, De Pontieu et al.| 
|2007[|Pereira et al.|2012| ). 
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